Dynamics of a hard sphere granular impurity 
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An impurity particle coupling to its host fluid via inelastic hard sphere collisions is considered. It is 
shown that the exact equation for its distribution function can be mapped onto that for an impurity 
with elastic collisions and an effective mass. Application of this result to the Enskog-Lorentz kinetic 
equation leads to several conclusions: 1) every solution in the elastic case is equivalent to a class 
of solutions in the granular case; 2) for an equilibrium host fluid the granular impurity approaches 
equilibrium at a different temperature, with a dominant diffusive mode at long times; 3) for a 
granular host fluid in its scaling state, the granular impurity approaches the corresponding scaling 
solution. 
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Granular fluids arc the subject of growing attention us- 
ing a variety of experimental, numerical, and theoretical 
tools 0, 0]- Among the many objectives is a theoret- 
ical characterization of important mechanisms based in 
the fundamentals of nonequilibrium statistical mechanics 
and kinetic theory 2]. In this context, an idealized model 
that captures the single feature of coUisional energy loss 
has been considered extensively. This is the system of N 
smooth, hard spheres with pairwise inelastic collisions. 
It has a remarkably rich number of properties that corre- 
late with those of real granular fluids. Of central interest 
are properties that differ from those of normal fluids, and 
the reasons behind these differences. 

Perhaps the simplest test bed for such questions is 
a study of the dynamics of a single impurity particle 
of mass Too in a- one-component granular fluid. The 
impurity-fluid particle interactions have a "restitution 
coefficient" measuring the degree of inelasticity (0 < 
ckQ ^ Ij with ao = 1 corresponding to elastic collisions). 
A complete description of impurity properties is given 
by the probability density for its position and velocity at 
each time, f (qo, vq, t). Its time dependence is coupled to 
the joint distribution for the impurity and a fluid parti- 
cle, /^^■' (xo,a;i,t), where Xi = (qi,Vi), through coUisions 
between the impurity and fluid particles. This coupling 
occurs through an exact equation 

(a* +V0 • Vo)^^(xo,<) = L[xQ-ao,mo\f^^\t)], (1) 

where L[a:o; ag, too|-] is a linear functional depending on 
the impurity position and velocity, and its restitution 
coeflicient and mass. Equation Q is the exact flrst 
Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hi- 
erarchy equation 0, with its detailed form given below. 
A primary result being reported here is the identity 

L[a;o; ag, toqI/^'^W] = L[a;o; «o = 1, mSI/^'^ (t)], (2) 



where the effective mass is 

toS = "T-o + (™o + (1 - ao) / (1 + ao) , (3) 

TO being the mass of a fluid particle. This identity es- 
tablishes an exact relationship of the coupling of the im- 
purity to its host fluid by inelastic and elastic collisions. 
Therefore, it has the potential to reveal details about the 
differences and similarities of granular and normal fluids. 

The main utility of Q is its basis for constructing a 
closed kinetic equation for F{xQ,t). In detail, the above 
functional L depends on J^^-* only for the colliding pair 
at contact, and on the pre-coUision hemisphere. For 
this configuration, a kinetic theory results by expressing 
Z*^^-* (xq, xi,t) as a linear functional of F{xo, t). A simple 
and accurate choice is the Enskog-Lorentz approxima- 
tion i 

/(2)(xo,xi,i) ^^^(xo,t)/(a:i,i)g(qo,qi|n(i)), (4) 

where f {xi,t) is the fluid probability density and 
9 (QO: ^.l W) is the local equilibrium pair correlation func- 
tion as a known functional of the nonequilibrium fluid 
density n (q, t). The correlation factor g accounts for the 
influence of the impurity particle on the joint impurity- 
fluid particle distribution. However, the one-particle fluid 
distribution / is not affected by the impurity in the ther- 
modynamic limit {N — > oo). 

Use of the approximation Q in the exact result 
gives the granular Enskog-Lorentz kinetic equation for 
F{xo, t). The utility of this equation is limited to low or 
moderate fluid density and impurity size small compared 
to the mean free path. However, within this context, its 
scope is quite broad, encompassing homogeneous and in- 
homogeneous states of the impurity and arbitrary host 
fluid states. The identity ((SJ implies that this kinetic 
equation is the same as that for an impurity with elastic 
collisions but with an effective mass. Several strong con- 
clusions follow immediately: 1) Since the effective mass is 
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a function of both the impurity mass and restitution co- 
efficient, every solution in the elastic case corresponds to 
a continuous family of granular solutions with different 
aQ,mQ. Equivalently, there is a scaling for the granu- 
lar solutions such that the dependence on aQ,mo occurs 
only through the single combination mg. 2) For the case 
of a host fluid at equilibrium, the granular particle ap- 
proaches equilibrium as well, but at a different temper- 
ature. In addition, the long time dynamics is hydrody- 
namic (diffusion) at long wavelengths. 3) For the case 
of a host granular fluid in its homogeneous cooling state 
(HCS), it is shown below that the impurity particle also 
approaches a HCS with a different temperature, but the 
same cooling rate. These conclusions hold for the com- 
plete range of inelasticity, < ao < 1- 

There are some limiting cases where some of the im- 
plications discussed here have been noted before. For 
a massive impurity particle the Enskog-Lorentz kinetic 
equation simplifies to a Fokker-Planck equation, and its 
exact mapping to the corresponding elastic case has been 
given 0. The equilibrium state of an impurity particle 
in an equilibrium host fluid has been discussed in 0. It 
has been noted recently ^ that the impurity mass can 
be taken as equal to the fluid particle mass by introduc- 
ing an effective restitution coefficient, while Piasecki et 
al. found (considering one-dimensional systems only) 
that the special case too/to — ao is equivalent to ao = 1 
with m^/m = 1. All of these special cases are subsumed 
by the general results discussed here. 

The system considered is a fluid of N hard spheres 
of mass m and diameter a, and an impurity particle of 
mass mo and diameter ctq. The impurity interacts with 
the fluid particles by inelastic hard sphere collisions; the 
nature of collisions between a pair of fluid particles is 
left unspecified at this point. The reduced distribution 
function for the impurity obeys the exact equation 
with the definition 

L[xo;ao,mo\f'-'^\t)] = a'^^'^ j dwida ■ ■ ^\ 

X ao V^^Hqo,v[,',qo -o^,v'/,t) 

-/'^^(qo,vo,qo + o^,vi,i) . (5) 

Here, = (tro 4- cr) /2 is an effective hard sphere diame- 
ter for an impurity-fluid pair, 'W = a a, d is the dimen- 
sionality of the system, Q is the Heaviside step function, 
Vol = Vo — vi is the relative velocity, and 
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(voi • 9)9 (6) 



are the restituting collision values that lead to {vo,vi} 
following a binary collision. It is easily verified that the 
total momentum of the colliding pair is conserved. Now 



make the change of variables Wi = Vi — (a^ — l) (voi ■ 
a)a in the first term of the integrand in ISJ, using dvi = 
aodwi. Vol • 9 = aowoi ■ 9, and Woi = vo — wi. After 
some straightforward algebra, the above scattering law 
becomes that for an elastic relationship between the pair 
(vo,Wi) and (v^',v'/), 

„ 2m _ 
Vq = Vo — (woi-cr)cr. 
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where mg (ao,mo) is defined by Q. Note that this does 
not imply equivalence between inelastic and elastic scat- 
tering laws. However, only the integration over Vi is re- 
quired in Q and this is equivalent to an integration over 
the new variable Wi. Finally, by renaming the dummy 
integration variable Wi vi the desired result, Eq. Q 
above, is obtained. 

There are some caveats to note at this formally exact 
level. The joint distribution /^^^ cannot be chosen arbi- 
trarily. In fact, it obeys an independent second hierarchy 
equation for which the dependence on ao cannot be re- 
moved by this transformation. Consequently, while the 
functional on the right side of |(5J is that for an elastic im- 
purity, in general its argument is not. Hence there is not 
a complete mapping of the granular impurity problem to 
one with elastic collisions. However, in the approxima- 
tion to be considered next this mapping is complete. 

Use of the approximate functional in the first hier- 
archy equation gives the revised Enskog-Lorentz kinetic 
equation 0,0] 

{dt + vo ■ Vo) F{xo,t) = lE[xo;ao,mo\F{t), f{t)]. (8) 

The collision operator, Ie, is a functional of 
both the impurity and fluid particle distribution 
functions, F{xo,t) and f(xi,t), with the notation 
lE[xo;ao,mo\F{t), f{t)] = L[xo;ao,mo\F{t)f{t)g]. At 
low density g — > 1 and this becomes the Boltzmann- 
Lorentz-Bogoliubov equation. The identity Q gives 
directly the corresponding property of the colli- 
sion operator lE[xo;ao,mo\F{t), f(t)] = lE[xo;ao = 
l,mo\F{t), f{t)]. Consequently, this imphes the equality 
of solutions, for any given fluid state, F{xo,t; ao,mo) = 
F{xo,t]ao = Ijmg). Thus, for any chosen positive value 
of mg (ao, mo) there is an equivalence class of solutions 
for different pairs (ao, mo). Furthermore, the solution for 
each different class is that for an impurity with elastic col- 
lisions. These are the main observations of the present 
work for kinetic theory. It is emphasized that they apply 
for all solutions, including time dependent and spatially 
inhomogeneous states, and for arbitrary nonequilibrium 
fluid states as well [lO|. In the following, some of the 
consequences are noted for illustration. 

For simplicity, the rest of the discussion is limited to 
cases for which the host fluid is in a spatially homo- 
geneous state, /(qi,vi,t) — > f{vi,t). Then, evaluated 
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at contact, 5 (qo, QiI'T- (*)) becomes an overall constant 
function x('^) of the time independent global density 
n = N/V. Aside from this factor, the collision operator 
becomes the usual Boltzmann-Lorentz operator for an 
impurity of mass uiq, i.e., Ie[xo, Q^o = Ij "^ol-^(^)' /(Ol ~^ 
x(n)/B[vo|F(t), /(t)]. Since the equation is linear it is 
sufficient to consider a single Fourier mode and the ki- 
netic equation becomes 

dtFik, V, t) = -tk ■ vi?(k, V, t) + x{n)lB [v|F (t) , f{t)]. 

(9) 

This is the equation for an impurity with elastic colli- 
sions, all effects of the impurity inelasticity being con- 
tained in the effective mass ttiq. 

The condition for a steady state F{v,t) -Fe(v) is 
k = and /^[vlFe, /(t)] = 0. For an isolated system 
(e.g., unbounded, or periodic boundary conditions) the 
only solution is a Maxwellian for both the fluid and im- 
purity distributions, at the same temperature, 

F,iv) = im*/2nT*f^\'"'>"/^^o\ (10) 

with Tq* = T. The distribution for the impurity particle 
can be expressed in terms of its actual mass mp by intro- 
ducing the true impurity temperature Tq = itiqTq /niQ — 
r(f +q:o)/[2+(1 — ao)m/mo]. Thus, the impurity particle 
is at equilibrium (Maxwellian), just as the fluid particles, 
but at a lower temperature Tq < T. This agrees with pre- 
vious results for the stationary state of an impurity in 
an equilibrium fluid. 

Next, consider the fluid state to be in equilibrium, 
f — fe, and study the approach of the impurity to its 
equilibrium distribution from a general initial state. For 
this purpose write 

^^(r,v,i)=Fe(v)[l-H0(r,v,t)], (11) 

so the Enskog-Lorentz equation becomes (in Fourier 
space) 

dt4>{k, V, t) = (-^k • V + Xe) ^(k, V, t), (12) 

where I^X = x('^)-'^e~^^s[v|FeX, /e]. The spectrum of 
—ik • V -|- has been characterized rigorously in the fol- 
lowing sense • The spectrum of le is negative with an 
isolated point at zero (conservation of probability). The 
spectrum of — ik ■ v -I- Xe is analytic in k about k = 0, 
and its real part is negative. Hence the point at zero be- 
comes one at —D{k) and remains isolated from the rest 
of the spectrum for sufficiently small k. In that limit 
D{k) Dk^ and represents a diffusion mode. In sum- 
mary, an initial state decays rapidly to a dominant dif- 
fusion mode and subsequently approaches uniform equi- 
librium by relaxation of this hydrodynamic mode. This 
is the classical picture of "aging to hydrodynamics" , fol- 
lowed by approach to equilibrium. The significance of 



this well-established result is that it now applies as well 
to a granular impurity, through only a change in the ef- 
fective mass. 

If the host fluid itself is composed of granular particles 
with inelastic collisions then the homogeneous state is no 
longer equilibrium but rather the HCS 0, 0| 

/(v, t) ^ nc-^{t)h {v/c{t)) , c{t) = ^2T{t)/m, (13) 

where T{t) is related to the average kinetic energy in the 
usual way. Its time dependence is due to the continual 
coUisional energy loss, or "cooling". A similar scaling 
solution for the impurity particle is sought of the form 

F(v, t) ^ c-^{t)Fh iv/co{t)) , coit) = ^2T*{t)/m*. 

(14) 

Here (t) is the corresponding measure of the impurity 
particle's average kinetic energy. Its time dependence is 
determined from the Enskog-Lorentz equation with an 
explicit time dependence through the host fluid temper- 
ature. If the impurity and fluid velocities are scaled with 
co(t) and c(t), respectively, and a dimensionless time s is 
defined through ds — [co{t)/i]dt, with i = 1 / xin)nW^^^ , 
then the Enskog-Lorentz collision operator becomes a 
function of time only through TQ{t)/T{t). The scaling 
solution (|14|l exists, therefore, only if this ratio is con- 
stant, T^{t)/T{t) = 7. This does not imply that the 
temperatures are the same, only that the cooling rates 
for the fiuid and impurity particle are the same. 

The specific form for Ff^ (w/co(i)), which is in general 
different from that of fh{v/c{t)), is determined from the 
Enskog-Lorentz equation. In the above dimensionless 
units, with the cooling rate denoted by (q = —ds In Tq (t), 
it becomes the time independent stationary equation 

CFh^O, /: = -^Vv-v+X. (15) 

Here, XX = /b[v|X, /^J, except that the dimen- 
sionless relative velocity is vqi = vq — viy^m^T^T 
and the dimensionless restituting velocities are Vq ~ 
Vq - [2m/ {ml + m)] (voi • ^)^ and v'/ = vi + 
[2y^mjw7/ (ttIq -|- m)] (vqi • In this context, it 

might be expected that initial states for the impurity 
approach this stationary HCS solution, analogous to the 
approach to equilibrium described above. Let us prove 
that this is indeed the case. 

Paralleling the above analysis, define 

T(r, V, t) = c-Q%t)Fu{v/co{t)) [1 + 0(r/£, v/co(t), t)] , 

(16) 

so that for a single Fourier component one has 

(-^k•v+Z)0, £ = -^vVv-(Xl)+X, (17) 

with £ = F;;^CFh and T = Ff^^IFh- Define a Hilbert 
space with scalar product (a, b) = J dvFh {v) (v) h (v). 
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Then, the following identity applies 

Re (V-, ZV') = Re (V', li^) - ^ (V, V-^l) • (18) 
Next, note the representation 

((/?, JV) = y dvodvidSF,,(vo)/h(vi)e(voi • 5^)|voi ■ B\ 
X [ip+ (v^') - (^+ (vo)] ^ (vo) , (19) 

where use has been made of the consequence of elastic 
collisions dv^dvi = dvgdv", voi • a- — — Vq^ • S, and the 
fact that the inverse and direct scattering laws are the 
same. Applying this to the right side of p8|) gives the 
desired result 

Re(7/>, (-ik-v + Z)V') <0, (20) 

i.e., the generator for the dynamics is dissipative [l^. 
A similar analysis of ((^, £-0) shows that a constant is 
an eigenfunction of both C and its adjoint, with zero 
eigenvalue. The null space is therefore spanned by 1. By 
normalization the solution (j) in (|17|l is orthogonal to 1. 
Consequently, 

a,||0f = 2Re(0,a,,^) <0. (21) 

Since the inequality is strict now, \\4>\\ is positive, mono- 
tonically decreasing, and must go to zero. 

In summary, the connection 10) obtained here presents 
the possibility to explore new consequences for the ki- 
netic theory of granular impurity dynamics. This has 
been illustrated for the case of the Enskog-Lorentz ki- 
netic equation. The result p2|) for the equilibrium host 
fluid is significant since it establishes the fact that in- 
elasticity does not affect any qualitative features of the 
impurity dynamics in that case. In particular, the exis- 
tence of a diffusion mode, its dominance at long times, 
and approach to equilibrium still apply. In the more in- 
teresting case of a granular host fluid, the result (|21|l is 
significant in showing that the non-trivial scaling solution 
Fh is universal in the same sense as the Maxwellian for 
normal fluids. Further study of the inhomogeneous gen- 
erator — ik • V -|- £ is needed to answer such questions as: 
do spatial perturbations exhibit a diffusive mode?; if so, 
does it exist and dominate at long times even for strong 
dissipation? Formally, the diffusion mode can be identi- 
fied by assuming spectral isolation and analyticity in k 
for a perturbation expansion of the zero eigenvalue of C 
|l4| . However, analysis in the large mass limit suggests 
that the answers may be conditional. In that limit it has 
been shown that diffusion requires that the cooling 
rate of the host fluid not be too large compared to the 
impurity-fluid collision rate. More complete characteri- 
zation of the spectrum of the generators of dynamics for 



inhomogeneous states seems to be limited at present to 
idealized models [l^ll5| . 

Kinetic theory for granular fluids and its domain of 
validity remain topics of intense debate. In practice, lim- 
itations due to approximate implementations can become 
confused with those of the theory itself. In this context, 
precise results for restricted systems (e.g., an impurity 
particle) can provide instructive benchmarks for impor- 
tant conceptual issues. This is one of the objectives of 
the present analysis. 
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